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We study how the thermodynamic properties of the Triangular Plaquette Model (TPM) are in¬ 
fluenced by the addition of extra interactions. The thermodynamics of the original TPM is trivial, 
while its dynamics is glassy, as usual in Kinetically Constrained Models. As soon as we general¬ 
ize the model to include additional interactions, a thermodynamic phase transition appears in the 
system. The additional interactions we consider are either short ranged, forming a regular lattice 
in the plane, or long ranged of the small-world kind. In the case of long-range interactions we call 
the new model Random-Diluted TPM. We provide arguments that the model so modified should 
undergo a thermodynamic phase transition, and that in the long-range case this is a glass transition 
of the “Random First-Order” kind. Finally, we give support to our conjectures studying the flnite 
temperature phase diagram of the Random-Diluted TPM in the Bethe approximation. This corre¬ 
spond to the exact calculation on the random regular graph, where free-energy and conflgurational 
entropy can be computed by means of the cavity equations. 


I. INTRODUCTION 

There are two main points of view to understand the 
nature of the glass phase. On one hand there is the idea of 
“dynamic facilitation”, which emphasizes the role of frus¬ 
tration in the dynamics: the motion of the microscopic 
constituents of a glass-forming system becomes inhibited 
by their close neighbours when the system is cooled. The 
consequence is that the material remains stuck for an 
extremely long time in a certain amorphous configura¬ 
tion. On the other hand there is the landscape scenario, 
or “Random First-Order Transition” (RFOT) theorjE® 
according to which the formation of a glass is the reflex 
of the existence of metastable states whose multiplicity 
is strongly reduced as the temperature is lowered. In the 
dynamic facilitation scenario thermodynamics is deem- 
phasized and the interaction only plays a role in dynam¬ 
ics. This idea is at the basis of the description of glasses 
provided by Kinetically Constrained Models (KCMJpESl^ 
KCM are lattice models where the variables do not have 
energetic interactions, but are subject to dynamic con¬ 
straints. Both theories have reference models that repro¬ 
duce important aspects of glass phenomenology, and it is 
hard to decide which scenario is the appropriate one to 
describe real systems. The two descriptions are very sim¬ 
ilar at the mean-field level: recent numerical simulations 
have shown that both at the level of average behavioi^ 
and at the level of fluctuationJi^ KCM models on ran¬ 
dom graphJi^ follow a glass transition pattern predicted 
by Mode Coupling Theory. IrfSl it was also shown that 
a typical RFOT model, the XOR-SAT on random graph, 
can be mapped in a KCM. The XOR-SAT model com¬ 
bines salient features of both theories: it has a trivial 
high-temperature thermodynamics as in KCM, but with 
a finite temperature entropy crisis glass transition as in 
RFOT. It is an interesting question if this commonality 


of mechanisms observed in mean-field extends to flnite 
dimensions. 

To investigate this question, in the present work we 
analyze the thermodynamic properties of a modified 2D 
Triangular Plaquette Model (TPM^PSHUI The original 
TPM is an example of KCMs: it is a spin model whose 
thermodynamics is the one of a trivial paramagnet while 
at the same time the model displays dynamical glassy 
phenomena with a super-Arrhenius relaxation time. In¬ 
deed, the TPM is nothing but a realization of the XOR- 
SAT model in finite (two) dimensionJSl. Our attention 
was brought to the TPM iim)articular by the results of^^. 
IiP, and more recently irP^, it has been shown that the 
TPM, in presence of external fields, supports both dy¬ 
namic and thermodynamic phase trans itions to glassy 
arrested phases. Emphasis is put liJUUol on the fact that 
as soon as such external fields are switched off ergodic- 
ity is restored, so that the TPM is just marginally glassy. 
Looking at the results of o f^^ l ^^ Rrom another perspective, 
they provide an evidence that the trivial thermodynam¬ 
ics of the TPM can be dramatically altered by means of 
very small perturbations. Following this line, our pur¬ 
pose here is it to show that the triviality of the TPM 
thermodynamics is marginal and its physics is close to 
the one of the landscape scenario. We will show that, 
as soon as some new interactions are introduced in the 
TPM, its thermodynamics cannot be trivial anymore and 
a phase transition appear in the system. Our results sup¬ 
port the point of view that the dynamic facilitation and 
the landscape scenario should be regarded as complemen¬ 
tary rather than alternative. 

Specifically, the TPM is a lattice spin model where 
the spins sit on a two-dimensional triangular lattice en¬ 
dowed with “plaquette” interactions: each plaquette cor¬ 
responds to the product of the three spins placed at 
the corners of an upward triangular cell of the lattice. 
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The energy of the model is H = — J2{ijk) where 

each triplet (ijk) of indices is associated to a plaque- 
tte. The remarkable observation made iiP is that there 
is a one-to-one correspondence between spins and pla- 
quettes: a configuration of the system is well defined 
either assigning the values of the spins or 

the values of the plaquette variables {ra}a=i,...,Af, where 
'I'a — (^i(a)'^j{a)(^k{a)- In terms of the plaquette vari¬ 
ables the Hamiltonian is equivalent to that of a system 
of non-interacting spins in a field, for which the partition 
function can be trivially calculated: Z = 2^[cosh(/3)]'^. 
The absence of thermodynamic singularities, to gether 
with a critical slowing down at low temperature j^^ * ^^ , 
Trei exp(H/T^), are typical features of KCMs. If one 
performs local MC updates acting upon the spins but 
then looks at the resulting dynamics of the plaquette 
variables, the latter looks like a kinetically constrained 
dynamics. It happens that the annihilation of a “de¬ 
fect”, namely the flip of a plaquette “a” from = — 1 
to Tq = 1, is favoured only when there is at least one 
other defect, i.e. excited plaquette, connected to “a”. 
Two plaquettes are connected when they have a spin in 
common. Since the dynamics acts upon spins, and in the 
TPM the update of one spin always corresponds to the 
updated of three plaquettes , the Monte Carlo dynam¬ 
ics corresponds to the flipping of three plaquettes per 
time. Due to the odd number of plaquettes attached to 
a spin (Ti, in the TPM there is no spin flip with AE = 0, 
since E = — X)aG9i’^“- same reason one needs 

that at least two of the plaquettes attached to ai are ex¬ 
cited, namely they are both = —1, in order to have 
AE < 0 by flipping ai. This explains why the anni¬ 
hilation of a defect is favoured only in the vicinity of 
another defect. The transition rates for spin updates 
depend on temperature through the standard metropo¬ 
lis rule, namely each attempted update is accepted with 
probability p = min(I, The idea we want to test 

in this paper is that if we introduce some new interactions 
between variables (or if we remove some), an important 
parameter in the thermodynamics is the ratio a = M/N 
between the number M of plaquettes in the Hamiltonian 
and the number N of spins in the systems. While such a 
scenario is well established for plaquette models on ran¬ 
dom graph^, this has not been tested to our knowledge 
in finite dimensional geometries. Of course the ratio a 
can be changed in many different ways. Here we will fo¬ 
cus on two class of models with extra plaquettes: in a first 
class we choose the triplets of spin in the new interactions 
completely at random, in a second class the new interac¬ 
tions are taken on a regular sublattice of the triangular 
lattice, either a fraction of them chosen randomly or all 
the interactions of the sublattice. With the first choice we 
induce arbitrarily long-range interactions, so that the re¬ 
sulting model is a kind of small-word networlP^. In order 
to be general, in the model with long-range interactions 
we also take into account a dilution of the plaquettes of 
the original two-dimensional triangular lattice: that is 
why we call such a model the Random-Diluted Triangu¬ 


lar Plaquette Model. The paper is structured as follows. 
In the first part, Sec. |llj we present the two class of mod¬ 
ified TPMs just mentioned. For each class we discuss the 
behaviour of the high-temperature expansion and present 
the results of numerical simulations. In Sec. IHII we dis¬ 
cuss how to use the leaf-removal algorithnPIl, which is 
a method borrowed from the study of constrained opti¬ 
mization problems, to draw a tentative phase diagram at 
T = 0 of the Random-Diluted TPM in the {as, ul) plane: 
tti is the concentration of long-range plaquettes while 
Us is the concentration of short-range plaquettes in the 
two-dimensional lattice. Finally, in Sec. |IV| of the paper 
we present the phase diagram at finite temperature for 
the Random-Diluted TPM on the random regular graph, 
which we will refer to often in the paper also as the Bethe 
lattice geometry, where the glass transition temperature 
can be exactly calculated by solving the cavity equations. 


II. TRIANGULAR PLAQUETTE MODEL WITH 
ADDITIONAL INTERACTIONS: 
HIGH-TEMPERATURE EXPANSION AND 
NUMERICAL SIMULATIONS 


The simplest choice of additional interactions for a 
modified TPM with Hamiltonian iJ = iFxPM + I^extra is 
represented by new ferromagnetic plaquettes: ILextra = 
— J2ijk^i'^j'^k- In the Hamiltonian iFextra the only 
source of randomness is then represented by the choice of 
which spins participate in each of the new plaquettes, i.e. 
the choice of the triplets of indices ijk. The ground state 
of the so modified TPMs is the configuration where all 
spins attain the value ai = -I-1. The aim of this work is to 
provide an evidence that in some of these modified TPMs 
a glass transition takes place before the system has the 
time to relax to the ordered ground state. We already 
know from the literature that this is the case for the TPM 
on a random graphP^. The idea of introducing new inter¬ 
actions is motivated by the purpose to induce an entropic 
crisis in the system, which is the typical mechanism for 
the formation of a glass phase within the RFOT theory 
scenaricP^l. We have already mentioned that in the TPM 
the partition function is simply Z = 2'^[cosh(/3)]^. Let 
us assume the existence of a modified TPM such that 
the number of plaquettes M is different from the number 
of spins N, but the partition function is still the trivial 
one: Z = 2'^[cosh(/3)]'^. In this hypothetical TPM the 
entropy of the system would be 

s{/3) = log(2) -h ^ log(cosh(^)) - ^/3 tanh(/3) (I) 


which in the limit of zero temperature yields 


lim s(/3) = ( 1 - ) log(2) 


/3—>-oo 


N 


( 2 ) 


Eq. tells us that in the case when a = M/N > 1 
the entropy is negative at T = 0, so that an entropic 
crisis takes place at T > 0. Clearly, since the TPM and 
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any decoration of it are models with discrete variables, 
any expression yielding a negative entropy cannot be an 
exact one. The only possibility is that, in a TPM with 
additional plaquettes {M/N > 1), the expression of the 
entropy in Eq. Q is an approximated one, for instance it 
can be the one provided by high-temperature expansion. 
This expansion usually well describes the properties of 
the liquid phase: the finding of a negative entropy by 
its prolongation to low temperatures tells us that the de¬ 
scription of the system as in the liquid state becomes 
inconsistent and the glass comes into plaji^. Our pur¬ 
pose here is to show that, in some modified TPMs, the 
parameter which controls the entropy crisis in the high- 
temperature phase is the ratio a = MjN between the 
number of plaquettes and the number of spins. While 
it is well know for mean-field plaquette models that the 
glass transition is controlled by such a ratio a = 
we provide here some evidences that this should be the 
case also in finite dimensions. All the TPMs with addi¬ 
tional interactions we discuss here are characterized by 
an Hamiltonian of the kind 


H T H^extrai 


( 3 ) 


where two different contributions read off: ^fextra = 
— wfth tho fndox r running over the new 

additional plaquettes, and the standard Hamiltonian of 
the TPM model, i?TPM = — Sfii j with the in¬ 

dex s running over the plaquettes of the two-dimensional 
triangular lattice. The two parameters which character¬ 
ize the model are the concentrations of original TPM 
plaquettes and of new additional plaquettes: respectively 
as = Mg/N and aL = Ml/N. The models with addi¬ 
tional plaquettes connecting arbitrarily far apart spins 
and those with additional short-range plaquettes will 
be presented respectively in Sec |H A| and Sec. |HB[ In 
Sec HA and Sec. HB we study the high-temperature ex¬ 
pansion and present some numerical results for the sim¬ 
plest case in which is arbitrary but ag is fixed to 
tts = 1, that is when there is no dilution of the plaquettes 
in the original 2D triangular lattice. Since the analysis 
of subsections Sec IH Al and Sec IHBI shows that the best 
candidate to display a low temperature glass transition 
is the TPM model with additional long-range plaquettes 
(“small-word” network), this model will be studied in the 
rest of the paper (Sec. HI and Sec. IVI also taking into 
account values < 1. The model in which both ag and 
ttL take values in the interval [0,1] is the Random-Diluted 
TPM. 


A. Long-range additional plaqnettes: the 
Random-Dilnted TPM 


A generalized TPM with additional plaquettes has the 
following high-temperature expansion for the partition 


function: 


Z = 2^ [cosh(/3)]'^ ^1 -I- ^ C(m, M)[tanh(/3)]™^ , 

(4) 


where the sum on the right hand side of Eq. Q runs over 
hyperloopTciade of m plaquettes: C(m, Mj represents 
the number of hyperloops of m plaquettes that can be 
built in a system with a total number M of plaquettes. 
An hyperloop is definecP^ as a set of interactions such 
that each spin is appearing an even number of times. The 
precise expression of C{m, M) depends on the model. For 
instance, in the pure TPM, one has C{m,M = N) = 0 
and the high-temperature series yields exactly the parti¬ 
tion function. Other situations will be discussed in the 
following paragraphs. Let us focus here on the behaviour 
of the high-temperature expansion in the case when the 
Ml = a^N plaquettes of Hextra have spins drawn with 
uniform probability from the lattice, with the only con¬ 
straint that spins in the same plaquette are different. Let 
us indicate with TZ an instance of disorder, namely a par¬ 
ticular choice of the triplets of spins in the plaquettes of 
Hextra- It is also Convenient to introduce the average of 
a function f[cr], with cr denoting a configuration of the 
system cr = (cti, ..., ctat), with respect to the measure 
provided by the pure TPM: 


(/M)tpm = ^ E W (5) 

^TPM „ 


The partition sum for a given TZ reads then: 


= 2'^[coi 


3sh(/3)]^+“^ tanh(/3))^ 


Ml 


= 2^[cosh(/3)]^+“^ 1 + E (a)tpm tanh(/3) -k 


rGTZ 


Ml 


+ E ('rr-A))TPM[tanh(^)]2-p... 


( 6 ) 


where, to lighten the notation, we used the plaquette 
variables Tr = cci^aj^ak^ and is the number of addi¬ 
tional plaquettes in iLextra- From Eq. (§ we see that 
the product of m “long-range” plaquettes contributes 
an hyperloop in the high-temperature expansion if and 
only if in the original TPM the corresponding corre¬ 
lation function is finite in 

the thermodynamic limit. Let us now notice that in our 
modified TPM with long-range plaquettes the annealed 
and quenched averages are equivalent at high tempera¬ 
ture due to the fact that the partition function is self¬ 
averaging, i.e. Z = Z"^: the proof of self-averaging is 
presented in Appendix A. This equivalence of annealed 
and quenched averages in the high-temperature phase 
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is typical not only of the the Random Energy ModeP^, 
the simplest model with an entropy-driven glass transi¬ 
tion like the one we expect in our modified TPM, but 
also of all the p-spin models fully-connected or on a 
random (hyper-)graph. In all these situations the be¬ 
haviour of the thermodynamic potentials in the high- 
temperature phase is obtained from the annealed free- 
energy / = —{l3N)log{Z). We are therefore allowed to 
average over disorder the partition function: 


2 : = 


Af{n 


M 


'^2n = 


TZa 


Ml\ At(7?,i) 


(Ml) 

= 2^[cosh(^)]^+^^ 1 -f ('^r-)TPMtanh(/3)-b 


/Ml\ 


'' ' r,p—l 


( 7 ) 


In Eq. 0 the indices r, p run over the set of all possible 
choices of random triplets of spins, while Af{TZk) is the 
number of ways in which k random triplets can be chosen, 
which, in the limit N ^ k ^ 1, reads 


A7(7^fc) 



( 8 ) 


In each sum on the right-hand side of Eq. Q the number 
of choices for the plaquettes not appearing in the brack¬ 
ets ( has been canceled out with the corresponding 
factor in the normalization constant. The binomial coef¬ 
ficient in front of each summation symbol in Eq. rep¬ 
resents the number of ways to choose m plaquettes out 
of Ml. Taking then into account that the typical value 
(TrTpTq )lpm does not depend on the choice of the indices, 
the expression in Eq. Q simplihes to 


Z = 2^[cosh(/3)]^+^^ • 


Ml 

m—l 


Ml\ 
m I 


(Ti...T„)^PM[tanh(/3)]" 


( 9 ) 


It can be then proven, as shown in Appendix A, that 


lim 

N—>-oo 



, .Trr 


= 0 , 


( 10 ) 


so that the high-temperature expansion is trivial and we 
have 


Z = 2^[cosh(^)]^+"^. 


( 11 ) 


The consequence of Eq. (11) is that, since 


a = {N + Ml)/N > 1, according to the expres¬ 
sion in Eq. ([^, at T = 0 the entropy is negative, which 
in turn implies an entropic crisis at T > 0. As already 
mentioned, the finding of negative entropy is, for a 
model like ours with discrete variables, an artifact of the 


high-temperature expansion: a phase transition usually 
takes place preventing this to happen. Such a phase 
transition is not necessary a glass one: indeed we find 
from the numerical simulations discussed below that a 
first-order transition to an ordered state is taking place. 
As shown iiPl 

for a TPM on random graph, the presence 
of an ordered ground state does not spoil the possibility 
of a glass transition. 

Numerical Simulations In order to test the behaviour 
of the TPM in presence of extra plaquettes with ran¬ 
domly chosen spins we realized numerical simulations of 
the model. We studied a lattice with triangular cells 
and the shape of a rhombus with L = 128 plaquettes 
per side, periodic boundary conditions, and a concen¬ 
tration of extra random plaquettes ol = 0.1. In order 
to speed up the equilibration dynamics at low tempera¬ 
tures we used a rejection-free algorithnP^. The behaviour 
of the internal energy along an hysteresis cycle is rep¬ 
resented by data in Fig. The scenario is as follows: 
by cooling down the system from high temperatures one 
first finds a spinodal temperature where an ordered 
metastable state appears. By further cooling the sys¬ 
tems it is then found a melting temperature where 
the ordered state becomes stable. The melting tempera¬ 
ture has been determined from the free-energy obtained 
numerically via thermodynamic integration: /(/3i) = 

PoPi^fWo) + dl3e{l3). In panel a) of Fig. 

there are two data sets which correspond respectively to 
the paramagnetic phase, at the higher energy eP“®'(/3), 
and the ordered phase, at the lower energy e^®™(/3): 
these two phases are stable respectively above and be¬ 
low the melting temperature T^. By independently per¬ 
forming the thermodynamic integration over eP®'''®(/3) and 
eferro(/3) one gets respectively /p®'''®'(/ 3) and /^“'’'°(/3): 
is determined by the crossing of /p®'“^^(/ 3) and /^®‘'"'°(/?). 
Since for T > 1 data are well interpolated by the high- 
temperature expansion, we fix the integration constant 
/{/3o) to /(/3o) = -/3olog(2) - /3o(l -b Q!LUog[cosh(/3o)], 
with /3(C^ = 1.2. In the top panel of Fig. IR the melting 
temperature is represented by a vertical dotted line: 
it is interesting to note that the system can be cooled 
to a remarkable extent below Tm while remaining in the 
disordered liquid phase. Bringing T down further, the re¬ 
laxation time becomes so large that the system falls out of 
equilibrium. Our cooling protocol is represented by runs 
ranging from 10® to 10^ Monte Carlo sweeps for each 
temperature, and by temperature jumps of AT = 0.04. 
Using this protocol we were not able to detect any ten¬ 
dency of the system to relax to the ordered ground state. 
From simulations we learn therefore that a first-order 
transition is present at T^,, but the systems is highly 
stable in the supercooled liquid-phase, namely at tem¬ 
peratures T < Tjn- This finding suggests that the system 
avoids on the time scales we sampled the negative entropy 
obtained by extrapolating the high-temperature expan¬ 
sion just by forming a glass. Let us also note that the 
annealed high-temperature expansion in Eq. © has a 










5 


very good agreement with simulations: in Fig. 0it can be 
seen that within the whole range of temperatures where 
we could equilibrate the system no relevant departure of 
data from the annealed energy e(/3) = —{1 + aL) tanh(/3) 
can be detected. In Fig. we can also see that the tem¬ 
perature To of the entropy crisis, as can be estimated 
from the high-temperature expansion of the entropy in 
Eq. 0, lies below the temperature range where we can 
equilibrate the system. With respect to the ideal glass 
transition temperature Tr-, the temperature Tq can be 
regarded as a lower bound for the possible values of Tx- 
It is hard to prove or disprove the existence of an ideal 
glass transition for our modified TPM solely on the basis 
of numerical data: this is clear looking, for instance, at 
the inset of the top panel of Fig. where the equilibrium 
relaxation time Tj-ei is shown as a function of tempera¬ 
ture. Good fits of the data (see the inset of the top panel 
in Fig. 0 can be obtained either with a function diverg¬ 
ing at finite temperature, r^ei ~ exp(A/(T — Tx)), or 
with the super-Arrhenius law Tj-ei ~ exp(i3/T^) (where 
clearly A ^ B)-. this is an ambiguous situation that is 
well know for these kind of model^^ and in general for 
glasses. That is why, in order to have a better insight 
into the thermodynamics of the TPM with random pla- 


quettes, we propose in Sec. Ill a different approach based 
on tools and ideas borrowed from the study of constraint 
satisfaction problems. 



FIG. 2. (Color online) Triangular lattice for the generalized 
TPM with additional short-range plaquettes on a regular sub¬ 
lattice made of equilateral triangles of side 3. The original 
TPMs plaquettes and additional interactions are shown as 
light gray and dark gray triangles, respectively. 



FIG. 1. (Color online) Main: Energy hysteresis cycle for a 
TPM model with additional long-range plaquettes. We sim¬ 
ulated the model on a lattice with the shape of a rhombus 
(see for instance the representation of Fig. with L = 128 
plaquettes per side and a concentration of additional plaque¬ 
ttes ttL = 0.1. (Red) Circles: cooling (cooling rate from 
10® to 10^ MCS per AT — 0.04); (Green) Triangles: heat¬ 
ing; Continuous (black) line: high-temperature expansion, 
e(r) = —{aL + l)tanh(l/r). The melting temperature Tm 
is indicated by the dotted vertical line while the temperature 
To of the entropy crisis of the high-temperature expansion is 
the filled (black) square. Inset: Circles are relaxation time in 
the supercooled liquid phase as function of T, and lines are re¬ 
spectively the two parameters fit with aexp{h/T) (continuous 
red line) and the three parameters ht with cexp{d/{T— Tx))- 
The critical temperature obtained from the fit is Tx ~ 0.374. 


B. Short-range additional plaquettes 

In the previous section we studied the effect of long- 
range random interactions on the high-temperature ex¬ 
pansion of the triangular plaquette model. It is then 
interesting to investigate what happens in presence of 
short-range interactions. Our choice for the new kind 
of short-range interactions has been guided by the pur¬ 
pose to reduce as much as possible the corrections to the 
high-temperature expansion. Since in the TPM the cor¬ 
relation of any group of three spins placed at the corners 
of an equilateral triangle of side 3 is zero, we defined 
an additional sublattice formed by equilateral triangular 
cells of that kind. The additional plaquettes appearing 
in i?extra are then chosen as all the plaquettes of this 
regular sublattice. In this case we know that at least the 
first term of the high-temperature expansion vanishes. A 
representation of this odd sublattice is given in Fig. 
where the additional plaquettes correspond to the spins 
placed at the corners of equilateral triangles with side 3. 
Let us then rewrite the partition function as a series in 
tanh(/3): 


Z = 2^[cosh(/3)]^+^^^ 


Ml 


1 -b^(r^)TPM tanh(/?) -b 

r 


Ml 

+ Xl^'^'-'^p)TPM[tanh(/3)]2 -b ... 
r^p 


(12) 
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where the sums in Eq. (12) run over all the plaquettes 
of the odd sublattice. At this stage the expression in 
Eq. @ is exact: it is just a way of rewriting the par¬ 
tition function. The difference between Eq. (12) above 
and Eq. Q in the previous section is that in Eq. (12) 
each term of the kind 


Ml 

(Ai ■•■r^jTPM[tanh(^)]'=, (13) 


accounts for the hyperloops made with all the possible 
choices of k plaquettes from the additional sublattice, 
while in Eq. 0 only the hyperloops coming from a given 
random choice of the additional plaquettes, denoted with 
TZ, were taken into account. Due to the regular structure 
of both the TPM lattice and of the sublattice of addi¬ 
tional interactions the existence of finite hyperloops can 
be proved along the following lines. Let us consider a lat¬ 
tice, TPM plus regular sublattice, of finite size, with open 
boundary conditions. Think for instance that our system 
is the portion of the lattice in Fig. which includes the 
4 additional plaquettes there represented and the 36 cor¬ 
responding plaquettes of the TPM, for a total number 
of M = 40 plaquettes. With open boundary conditions, 
the number of spins in this system is N = 49. In this 
case nothing forbids the partition function to be exactly 
Z = 2^[cosh(/3)]^, since at zero temperature the corre¬ 
sponding entropy is positive: s{T = 0) = 1 — 40/49 > 0. 
In this case there is no constraint implying the presence 
of hyperloops: it might be that for a lattice of this size 
there is no finite hyperloop. If we think in general to a 
lattice with the shape of a rhombus, with the same num¬ 
ber L of TPM plaquettes on the horizontal and oblique 
side (consider Fig. to have an idea), and with open 
boundary conditions, its partition function in absence of 
hyperloops is 

Z = 2(^+1)" [cosh(/3)]^"+(^/3)" ^ 


Due to the open boundary conditions, if the number of 
plaquettes on each row is L the number of spins on the 
same row is L -|- 1. On the other hand there are L/3 
rows of additional plaquettes, each with L/3 plaquettes. 
This explains why the total number of spins is (L -I- 1)^ 
and why the total number of plaquettes is + (L/3)^. 
For simplicity, let us consider just the case when L = 
3m, where the integer number m indicates therefore the 
number of plaquettes from the additional sublattice in a 
row: the zero-temperature entropy per degree of freedom 
corresponding to the free energy in Eq. (14) is 


s(T = 0) = 1 - 


lOm^ 


(3m -I- 1)^ 


(15) 


From Eq. (15) we find the maximum value of the param¬ 


eter m for which is possible to not have hyperloops in the 
partition function: m* = 6. For m > m* the expression 
of s(T = 0) in equation Eq. © is negative: this actu¬ 
ally proves (the system has discrete variables) that for 


any lattice of size m > m* there must be hyperloops in 
the system. This way of reasoning allows one to prove not 
only that finite hyperloops are present, but also to fix an 
upper bound for the size of the smallest hyperloop. Since 
in a lattice with m* -1-1 = 7 there is necessarily an hyper¬ 
loop, then the smallest hyperloop cannot have more than 
(m* -|-1)^ = 49 plaquettes taken from the additional sub¬ 
lattice. Let us now explain why in this case there are non 
trivial contributions to the high-temperature expansion. 
In order to be general, we assume the smallest hyper¬ 
loops is made of k additional plaquettes of the sublattice 
and by m plaquettes of the original TPM lattice. For 
every correlation function of k' plaquettes such that k' 
is a multiple of k, and the plaquettes are all taken from 
the sublattice of additional ones, one must take into ac¬ 
count disconnected hyperloops. The leading contribution 
of this kind of terms to the high-temperature expansion 
is 


Z = 2^[cosh(^)]^+^^ 


1 -I- ML[tanh(/3)] 


fc+m 


+ Ml{Ml - fc)[tanh(/3)]2'=+2™ + 

+ Ml{Ml - k){ML - 2A:)[tanh(/l)]3'=+3'" . 


(16) 


where = oilN. An explicit expression for the com¬ 
binatorial prefactor C(m,M) appearing in E q. (E] ) can 
be provided for the terms appearing in Eq. ( 16|). Let 
us indicate with C(mTPM + Wgxtraj Af) the combinatorial 
factor which accounts for the multiplicity of an hyperloop 
with mxpM plaquettes from the original TPM model and 
Wextra additional plaquettes. Then, according to the ex¬ 
pression in Eq. (16) and taking n a positive integer we 
can write 


C(A:n + mextra,M)~M". (17) 

The series in Eq.([I^ can be summed leading to 

Z ~ 2^[cosh(,5)]^+^^e“^^[‘“‘'(^)l'“""‘, (18) 

which represents the contribution to the high- 
temperature expansion provided by fc—plaquettes 
connected hyperloops. We can conclude that in the 
TPM the addition of short-range additional plaquettes 
always introduce corrections to the high-temperature 
series: it is then hard to say whether or not this 
expansion, which well reproduces the properties of the 
system in the liquid phase, has positive entropy at all 
temperatures. What numerical data (see below) show 
is that even short-range additional interactions induce 
a first-order transition with a rather robust supercooled 
liquid phase surviving below the melting temperature. 
In order to characterize the hyperloop corrections to 
the liquids phase, we looked to the smallest hyperloop 
for the model defined at the beginning of this section. 
In order to do this we used a simulated annealing 
method, explained in detail in Appendix] A 2| According 
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FIG. 3. (Color online) Smallest hyperloop of the high- 
temperature expansion for the modified TPM with extra 
short-range plaquettes discussed in Sec. |II B[ In the picture 
can be counted 54 small plaquettes of the original TPM lat¬ 
tice, thick (red) triangles, and 10 plaquettes of the auxiliary 
sublattice, thin (black) triangles, which has as elementary 
cells equilateral triangles of side 3. 


to simulated annealing, the smallest hyperloop of the 
high-temperature series is the one shown in Fig. 
which is made of 54 plaquettes from i/xPM and 10 
plaquettes from Tfextra- Even if such hyper loop is “big”, 
in the sense that it provides a correction of the order 
[tanh(,5)]®^ to the high-temperature expansion of the 
free-energy, it contains a number of plaquettes from 
the additional sublattice (10 plaquettes) well below 
the threshold of {m* -1-1)^ = 49 plaquettes (from the 
additional sublattice) discussed above in this paragraph. 


The next question can then be: what about if we 
add to the TPM Hamiltonian not all the plaquettes 
of the sublattice but just a fraction of them randomly 
chosen? That is, what about if the sums in Eq.(12) runs 
over the plaquettes of a special instance of the disorder, 
which corresponds to have in iJextra just a fraction c 
of the plaquette of the regular sublattice? The answer 
is that hyperloops will still be there, just in a smaller 
amount compared to having all the plaquettes of the 
regular sublattice. Let us consider for instance the 
hyperloop of Fig. made of 10 close-by plaquettes of 
the regular sublattice and 54 plaquettes of the TPM. For 
a random choice of the additional plaquettes such that 
a fraction c of the total number present in the sublattice 
is ta ken , in a high-temperature expansion like the one in 
Eq. (16) a contribution as ML[tanh(/3)]®^ is replaced by 
cMx[tanh(/?)]®^, and correspondingly the correction to 
the partition function becomes exp(caxA^[tanh(/3)]®'^). 
The difference between adding randomly chosen long- 
range and short-range plaquettes is at this point clear: 
while in the former case it can be proved that the 
high-temperature expansion is trivial and there is an 
entropic crisis, in the latter case there are corrections 


coming from finite hyperloops, and these corrections 
are most likely preventing any entropic crisis. For this 
case we present the numerical evidence that, even with 
no clue on the presence of an entropic crisis, the idea 
to add interactions to the TPM is a good strategy to 
induce a non-trivial thermodynamics characterized by 
the presence of a supercooled liquid phase. 


Numerical Simulations In order to provide the nu¬ 
merical evidence that even short-range interactions pro¬ 
duce a non-trivial thermodynamics we considered the 
case of a random choice of the plaquettes on the 
regular sublattice and we found convenient also to 
put weaker interactions on these additional plaque¬ 
ttes. Actually we considered the Hamiltonian iJextra = 
with Ji = 1/5 and with the sum 
running on a finite fraction of the sublattice plaquettes. 
On one hand the random choice of a subset of the pla¬ 
quettes on the regular sublattice decreases the number 
of finite hyperloops in the high-temperature expansion, 
on the other hand a reduced strength of interactions on 
these plaquettes decreases the weight of these hyperloops. 
Monte Carlo simulations done with Ji = 1 showed that 
also with short-range plaquettes there is a melting tem¬ 
perature Tra, but differently from the case of long-range 
plaquettes (with a random choice of the spins) the system 
decays to the ordered ground state as soon as T < Tm- 
Data are not reported here, but we found that with the 
same annealing protocol used for long-range plaquettes, 
namely temperature jumps of AT = 0.04 each 10® or 
10^ MC steps (depending on the temperature), the sys¬ 
tem reaches the equilibrium ground states for all temper¬ 
atures T < Tm with M]^ = O.IN additional short-range 
plaquettes. Before comment the numerical data obtained 
with Ji = 1/5, let us notice that when the additional 
plaquettes are taken from a regular sublattice, even for 
T > Tm the high-temperature expansion yields a small 
but finite magnetization. This can be clearly seen looking 
at the high-temperature expansion of the magnetization: 



[cosh(,d)]*^^ 

NZ 


_ ^ ^ p-/3(ffTPM[o-]-|-ffextra[CT]) _ 

i—1 cr 

~ N Ml 

^ ^((Ti tanh( Ji/3) -|- . . . 

_i—l r=l 


N 


tanh®(/3) tanh( Ji,5) -I-..., 


(19) 


which, as can be seen in Fig. nicely compares with 
the results of numerical simulations. As is clear from 
Eq. (19), the first effect of taking Ji < 1 is to reduce 
the high-temperature magnetization in the vicinity of Tm 
making the supercooled liquid phase more stable. Let 
us briefly explain how the equalities in Eq. (19) are ob¬ 
tained. From the first to the second line of Eq. (191 we 
just wrote down explicitly the high-temperature expan¬ 
sion. The last identity of Eq. (19) is clear if one recalls 
that each plaquette of iLextra is represented by a triplets 
































FIG. 4. (Color online) Main: Energy hysteresis cycle for a 
TPM model with additional short-range plaquettes. We sim¬ 
ulated the model on a lattice with the shape of a rhombus 
(see for instance the representation of Fig. with L = 128 
plaquettes per side and a concentration of additional plaque¬ 
ttes aL = 0.1: the coupling coefficient of the additional pla¬ 
quettes is Ji = 1/5. (Red) Circles: cooling (cooling rate 
from 10® to 10^ MCS per AT = 0.04); (Green) Triangles: 
heating; Continuous (black) line: high-temperature expan¬ 
sion, e(r) = —(ol tanh(Ji/r) -|-tanh(l/T)). Inset: Magne¬ 
tization hysteresis cycle for the same model with the same 
cooling/heating protocol. (Red) Circles: cooling; (Green) 
Triangles: heating; continuous line: m{T) from the high- 
temperature expansion: m(T) = Q:i[tanh(l/T)]® tanh(Ji/r), 
with ul = 0.1 and Ji = 1/5. 


of spins vertices of an equilateral trian¬ 

gle of side 3. It can be easily realized that those spins 
form an hyperfieltP^ when multiplied with the plaque¬ 
ttes of i/xPM which are enclosed within the perimeter 
of the same triangle. Borrowing the terminology oP^we 
call hyperfield a set of interactions such that all spins 
but one, ai, appear an even number of times. Since 
the plaquettes of ^/^extra are all represented by triplets 
of spins placed at the corners of equilater triangles of 
side 3, each of them contributes an hyperfield ai in the 
high-temperature expansion, and this hyperfield becomes 
an hyperloop of weight [tanh(,5)]® tanh( Ji/3) when mul¬ 
tiplied for the corresponding spin Ui appearing in the 
definition of the magnetization. Then, by taking also the 
approximation Z ~ 2^[cosh(/3)]^+'^^ (which is exact in 
the case of long-range plaquettes, see Sec. II A), we have 
the result of the last line in Eq. (19). 

Our numerical study showed that for the smaller value 
of the coupling constant Ji = 1/5, a number = O.IA^ 
of additional plaquettes and the same annealing proto¬ 
col already mentioned above in this section, at all the 
temperatures studied the system does not decay to the 
ordered ground state, see Fig. In the inset of Fig. 
is shown the behaviour of the magnetization, which in 
the paramagnetic phase is well reproduced by the high- 


temperature expansion. The model with randomly cho¬ 
sen plaquettes on a regular sublattice seems therefore a 
good one to reproduce the standard scenario of realistic 
glass-formers: a first-order phase transition to an ordered 
ground state plus a long-lived supercooled liquid phase. 
A system where a similar behaviour i s foun d is the Cou¬ 
pled Two Level System (CTLS) modeP^ESl, xhe CTLS is 
a non-disordered plaquette model which presents, as our 
TPM with additional plaquettes, a first-order transition 
to an ordered ground state and a metastable supercooled 
liquid phase. To summarize, we can say that even when 
is not possible to argue about the existence of a ther¬ 
modynamic glass transition, the addition of short-range 
plaquettes to the TPM induce a first-order phase tran¬ 
sition and the formation of a robust supercooled liquid 
phase. We can therefore argue that even with short-range 
interactions the presence of a non trivial thermodynamics 
is controlled by the ratio a = M/N between the num¬ 
ber of plaquettes and the number of spins. In plaquette 
models with short-range interactions the glass-forming 
ability is then a matter of competition between two time 
scales: the time scale to nucleate the the crystal and the 
relaxation time of the supercooled liquid. For a detailed 
discussion on how a stable supercooled liquid phase can 
be obtained by appropriately tu ning the cooling rate pro¬ 
cedure let us refer the reader tcPSESESl 


III. RANDOM-DILUTED TPM: PHASE 
DIAGRAM AT T = 0 


The analysis of the present section is based on the deep 
connection between the thermodynamics of plaquette 
models and the properties of solutions of the XOR-SAT 
problem. The latter is a constraint satisfaction problem 
which has been very successfully described within the 
landscape scenario of the ideal glass transitioiJ^. The re¬ 
lation between the TPM and the XOR-SAT comes from 
the fact that the ground states of a TPM can be obtained 
as the solutions of a XOR-SAT. By assuming the change 
of variables ai = (—1)”% a ground state of the TPM can 
be always written as the solution of the system of linear 
equations: 


^ n J =0, V p (20) 

/ mod2 


where with the notation i € dp we indicate the spins in 
the plaquette p. The ground states of the TPM are rep¬ 
resented by all the solutions of the system in Eq. (20), 
which is made of M linear equations, the plaquettes, in 
N variables, the spins. The topology of the interaction 
network is specified by the equations of the linear sys¬ 
tem. For instance, in the triangular 2D lattice of the 
TPM each spin Ui appears in three equations and in 
each equation it is coupled to the two other spins be¬ 
longing to the same plaquette. When the network formed 
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by spins and plaquettes corresponds to a random hyper¬ 
graph (which is a locally tree-like network with loops of 
order log(A^)) it is known that the properties of the so¬ 
lutions of the XOR-SAT are fully determined by the pa¬ 
rameter a = For a TPM on the random graph, 

even with non-disordered interaction^^, two transitions 
are found varying a, respectively at ad and Ounsat- When 
a < ad the system of equations in Eq. (201 is solvable and 
happens that few spins flips are sufficient to go from one 
solution to another: the set of solutions forms a unique 
cluster. On the contrary, when ad < a < Ounsat the 
set of solutions splits in clusters separated by extensive 
barriers: it is necessary to flip an extensive number of 
spins to pass from one cluster to the other. Finally, 
when a > Ounsat the system in Eq. (20) is no more 
solvable: at Ounsat takes place the SAT/UNSAT tran¬ 
sition. With respect to any random update algorithm 
designed to move across the phase space of variables 
the SAT/UNSAT transition represents an 
ideal glass transition. The non-satisfiable phase of the 
XOR-SAT coincide with the glass phase of the related 
TPM: as a consequence, when a > Ounsat the TPM (on 
random graph) is in the glass phase at T = 0, which 
means that it has an ideal glass transition at T > 0. 


A. Leaf-removal algorithm and T — 0 phase 
diagram 

The “leaf-removal” algorithm is a decimation scheme 
used to study the satisfiability of the XOR-SAT model 
on random regular graph^^. We discuss here how this 
algorithm can be used to investigate the properties of 
our Random-Diluted TPM. It is called “leaf” every spin 
which appears in only one plaquette (i.e. every variable 
appearing in a single equation of the XOR-SAT prob¬ 
lem): “leaf-removal” is a prescription to remove itera¬ 
tively from the graph all the spins which are (or become) 
leaves. The procedure is iterative because after the re¬ 
moval of each leaf new leaves may appear in the system. 
When a < ad leaf-removal is able to remove all the spins 
from the graph, while for a > ad the algorithm stops 
leaving the so called “core”, i.e. a set of spins among 
which no one is a leaf. The clustering of solutions of the 
XOR-SAT at ad corresponds to the formation of the core. 
The SAT/UNSAT transition takes place when the num¬ 
ber of equations (plaquettes) in the core, Me, becomes 
larger than the number of variables left on it, Nc- The 
critical value Ounsat can be determined by studying the 
ratio Mc/Nc = Oc in the core, and corresponds to Oc = 1. 
On a random regular graph the dependence of Me and 
Ne on a can be determined analytically in the thermo¬ 
dynamic limilP2. The leaf-removal algorithm can be used 
to study the formation of the core and the behaviour of 
Oc on the core for a XOR-SAT on every kind of topology. 
Nevertheless, only for random graphs it is proven that 
the formation of the core corresponds to the clustering 
transition and the value Oc = 1 to the SAT/UNSAT tran- 


sitiorPil. On finite dimensional topologies a proof of this 
correspondence is still lacking. This notwithstanding, we 
studied numerically the action of the leaf-removal on our 
Random-Diluted TPM, proposing a “tentative” phase di¬ 
agram. We compare this numerical phase diagram with 
that obtained by analytically solving leaf-removal for the 
representation of our Random-Diluted TPM on the ran¬ 
dom graph: in this case the leaf-removal analysis yields 
exactly the thermodynamic properties of the systerrf^. 
Let us stress that by running the leaf-removal algorithm 
in a finite-dimensional geometry, we may find both a crit¬ 
ical value a* for the formation of the core and a critical 
value a** where Me/Nc = 1, but there is not proof that 
they correspond respectively to the clustering and UN¬ 
SAT transition. This is the reason why the phase dia¬ 
gram of Fig.j^is just tentative. The conjecture that even 
in finite dimensions at a** really takes place the UNSAT 
transition, namely a** = Ounsat, is supported within our 
analysis just by the agreement we find between the nu¬ 
merical finite-dimensional and exact mean-field predic¬ 
tions on the phase diagram. 

Until now we have discussed about a single parameter 
a, but for the Random-Diluted TPM and its representa¬ 
tion on the random graph we need two: as , which repre¬ 
sents the concentrations of plaquettes in the 2D lattice, 
and a^, which represents the concentration of long-range 
plaquettes. The Random-Diluted TPM corresponds in 
practice to a random-graph structure, parametrized by 
tti, built on the top of a diluted two-dimensional net¬ 
work, characterized by a dilution parameter as- The 
probability that a spin is attached to I short-range pla¬ 
quettes depends on as and is: 

p,{as)={^a\(\-asf-\ (21) 

while the probability that a spin is attached to i long- 
range plaquettes depends on a^ and reads: 

= ( 22 ) 

The probability that a spin is attached to overall I pla¬ 
quettes is 

^0 

nn = A/'”^ ^ Pr[pts)vi-r{piL), (23) 

r—0 

with 

oo Iq 

A/" = Pr{as)Pi-r{aL) = 

i—Q r—0 

3 3 3 k 

= 1-EE Pri^s^Pk — ri^k) E EE Pr{o^s)Pk—r{^L'), 

r—0 k—r k—0 r—0 

(24) 

where £o = minjl', 3}. When the Random-Diluted TPM 
is represented on the random graph the small-word topol¬ 
ogy, i.e. the two-dimensional lattice plus few long-range 
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FIG. 5. (Color online) Panel a)-. “T = 0” phase diagram of 
the Random-Diluted TPM (see Eq. (23l). Symbols represent 
numerical data obtained running the leaf-removal algorithm 
on finite dimensional geometries, lines the results of analytical 
calculation on the Bethe lattice. Circles (blue) and the dotted 
line represents the critical line for the formation of the core, 
squares (red) and the continuous line are the critical line for 
the SAT/UNSAT transition. Panel b)-. Diamonds, numerical 
estimate in finite dimensions of the critical value ul for the 
SAT/UNSAT transition as function of \/iV , with N the size 
of the system, at fixed Os = 1; dotted line, linear fit of data. 
Panel c)\ zoom of the SAT-UNSAT transition line close to 
the point (a^, = !,«_£, = 0). Triangles, analytic result on 
the random graph; Squares, numerical data for a system with 
N = 10® spins; Dotted line, linear fit of analytic prediction 


connections, is lost. The only ingredient of the original 
model which is kept is the presence of two kinds of pla- 
quettes, each characterized by a different probability for 
the connectivity with the spins. For the representation 
on the random graph the adjectives “short-range” and 
“long-range” are therefore only conventional: the former 
denotes the plaquettes attached to spins with the proba¬ 
bility of Eq. (211, the latter plaquettes attached to spins 
with the probability of Eq. (22). 

Fig.0 shows the phase diagram of the Random-Diluted 
TPM obtained by running the leaf-removal algorithm in 
finite dimensions and that obtained by solving the cor¬ 
responding equations (see Appendix E in the mean- 
field (random graph) approximation. In both cases we 
determine a “critical” line (ttg”®, a™’'®) for the forma¬ 
tion of the core (clustering transition on the random 
graph) and a critical line (a™®®'*, from the con¬ 

dition Mc/Nc = 1 (unsat transition on random graph). 
Eirst of all, let us note the agreement between numeri¬ 
cal and analytical predictions on the location of the line 
(Q,unsat,^unsat) Since with Os = 0 the Ran¬ 

dom Diluted-TPM is perfectly equivalent to the XOR- 
SAT*^, for Os = 0 we recover the random graph re¬ 
sult = 0.91^^. By looking at the left part of 

the phase diagram in Pig. we are indeed not surprised 


that for Oi ^ 1 and Os <C 1 the analytic predictions 
on the random graph are in agreement with the numer¬ 
ical analysis in finite dimensions: the Random-Diluted 
TPM is almost a random graph for these values of the 
parameters. What is more surprising is to find an agree¬ 
ment between the numerical and the analytical estimate 
of (a™®®'*, in the bottom right part of the phase 

diagram in Fig. where <C 1 and Og ~ 1. In 
this region the Random-Diluted TPM is almost a two- 
dimensional model, while the analytical predictions are 
given for a random graph geometry. It really looks like 
that as soon as few long-range interactions are added 
to the TPM its physics becomes immediately well repre¬ 
sented by the random-graph. If the true thermodynamics 
of the Random-Diluted TPM could be predicted on the 
basis of the analytic results of Fig. we would say that 
as soon as any finite concentration of random pla¬ 
quettes is added to the TPM the model develops a finite- 
temperature glass transition. We do not still have a proof 
of that, therefore by now this is just a conjecture sup¬ 
ported by the agreement between numerics and analyt¬ 
ics on the location of the line. The small 

panels of Fig. illustrate the finite-size scaling analysis 
needed to assess the agreement of analytical and numer¬ 
ical predictions on the behaviour of the (a™®®'*,a™®®'*) 
transition line close to the point (a^ = 1 , 0 ^ = 0). 
Concerning the analytical and numerical data on the lo¬ 
cation of the (a®”®, a™''®) line in the bottom right part 
of the phase diagram of Fig. [^we find a certain disagree¬ 
ment: two main comments on this are in order. First, it 
is well known that the clustering transition which takes 
place at ad on the random graph, and which corresponds 
to dynamical ergodicity breaking, is a purely mean-field 
phenomenon which turns into a crossover in finite dimen¬ 
sions. From this point of view we are not concerned about 
the disagreement found between the analytical (mean- 
field) and numerical (finite-dimensional) results on the 
position of the line (a™®®, On the other hand, the 

numerical results on the location of the line (a™®®, a™®®) 
are per se interesting and represent a useful source of in¬ 
formation on the model, as will be discussed in the next 
section Sec IIIIBl 


B. Leaf-removal and exact calculation of Z 

The numerical estimate of the position of the line 
(a®™®, q;^”®) in the phase diagram of Fig.j^is quite inter¬ 
esting even in the case of a finite dimensional geometry. 
In this case the line (a®®''®, a™®®) is the upper boundary 
of the region where the partition function of the Random- 
Diluted TPM can be exactly calculated and has the ex¬ 
pression 

Z = 2^[cosh(/3)]^'’+“® = (2 [cosh(/3)]“'’+“^)^ . (25) 

The exact explanation of why leaf-removal allows us to 
check whether or not the partition function can be ex- 
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FIG. 6. (Color online) Panel a): Circles are the critical val¬ 
ues (as, a_t) for the formation of the core extrapolated in the 
thermodynamic limit from numerical data at finite N. The 
continuous line represent a quadratic fit of data y = ax^, with 

Panel b): 


IIIB 


prefactor a = 0.877, see also Eq. (261 in Sec. 
critical values aL{N) for core formation as unction of the 
size N of the system for different values of (different sym¬ 
bols), from top to bottom: Os =0.93, 0.95, 0.97, 0.98, 0.99. 
Panel c): concentration of spins in the core as function of ql 
for a fixed value of as = 0.8, with N = 1024^. 


actly summed, yielding the expression in Eq. (251, can be 
found in Appendix]^ The behaviour of the (a™™, 
line (numerical data) close to the point (os = 1 ,q;l = 0) 
in the phase diagram of Fig. which correspond also 
to the data in the main panel of Fig. shows that, 
among all the two-dimensional models belonging to the 
line (as,0), the only one such that any concentration 
ol > 0 of long-range plaquettes is “critical” is the orig¬ 
inal TPM. As “critical” we mean that the possibility 
to exactly compute Z according to the expression in 
Eq. (251 is spoiled as soon as any finite concentration 


of long-range plaquettes, ul > 0, is introduced. In order 
to make ourselves really sure about that we need to know 
the behaviour of the system in the thermodynamic limit. 
In the top right panel of Fig. is presented the study of 
finite size effects, while in the main panel of Fig. ap¬ 
pears the resulting estimate for some points of the line 
(a®”®, a™’'®). These point can be well interpolated with 
a parabola: 


a, 


-f C>(Af-i/2) ^ (1 _ 


core\2 


(26) 


The parabolic ht of the data in the main panel of Fig. 
represents the main evidence that the only model on 
the line {as,aL = 0) such that the addition of extra 
plaquettes is critical is the TPM model. To conclude 
this section let us notice that also in an almost finite¬ 
dimensional geometry, as = 0.8 and ul ~ 1, we find 
that the formation of the core is a discontinuous process, 
see panel c) of Fig. as is usually found for random 


geometries. This result supports the view that the 
physical properties of the Random-Diluted TPM can be 
well represented even on a random graph. Our finding 
that the formation of the core happens discontinuously 
upon changing ol also in finite dimensions is intriguing: 
it is not the first time that a similar discontinuous 
transition has bee n obse rved in finite dimensiond^SHaej 
The spiral model oP^^^, which is a KCM, is the exam¬ 
ple of a finite-dimensional system where a dynamical 
transition really takes place and is due to the formation 
in the system of an infinite compact cluster of “frozen” 
(i.e. not allowed to move due to the kinetic constraint) 
spins. A cluster of spins which cannot flip due to a 
kinetic constraint is not exactly reconducible to the 
leaf-removal “core” of our Random-Diluted TPM. Yet, 
if the mean-held scenario is predictive also for the 
behaviour in hnite-dimensions, the formation of the core 
should take place when, in order to move in phase-space, 
we need to hip an extensive amount of spins, which we 
may very roughly think about as a “frozen cluster”. 


IV. RANDOM-DILUTED TPM ON RANDOM 
GRAPH: PHASE DIAGRAM AT T > 0. 


In this section we present results on the hnite- 
temperature phase diagram of the Random Diluted-TPM 
model on the random regular graph (Bethe lattice), 
where the temperature Tk of the ideal glass transition 
can be computed exactly. 

According to the presence of both long and short-range 
plaquettes in the Random-Diluted TPM, the cavity equa¬ 
tions (Appendix 0 for its representation on the ran¬ 
dom graph are written by means of two diherent cav¬ 
ity helds, as is usually done for small-word network^^. 
The held Ua^i determines the probability distribution 
p(CTi) ~ all iiig plaquettes attached to ai 

but a are removed, and a is a long-range plaquette. 
The held Va^i determines the probability distribution 
p(o’i) ~ when all the plaquettes attached to ai 

but a are removed, and a is a short-range plaquette. The 
cavity equations, written and discussed in Appendix 
allow us to hnd the equilibrium values of the helds u and 
ri, from which all the thermodynamic potentials can be 
calculated (formulas are in Appendix]^. For hxed val¬ 
ues of Qfs and ahi the glass transition temperature Tk 
is obtained as the temperature where the conhgurational 
entropy E vanishes. In Fig. is represented the phase di¬ 
agram of the Random-Diluted TPM in the plane (a^, T) 
for two values of the short-range plaquettes concentra¬ 
tion: Os = 1 and as = 0.96. While for Og = 1 the 
glass transition temperature vanishes when also the con¬ 
centration Q!l of long-range plaquettes vanishes, when 
as = 0.96 we hnd that Tk vanishes at a hnite value 
of the additional plaquettes concentration, a™'” > 0: 
when a^ belongs to the interval [0, a™™] the system is 
liquid at all the temperatures. By dehnition, the hnite- 
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FIG. 7. (Color online) Phase diagram of the Random Diluted- 
TPM model on the Bethe lattice in the (a_L,T) plane for two 
different values of dilution Circles (red), Os = 1; Squares 
(blue), Us = 0.96; Continuous lines; fits of the data with the 
function aL{T) = Ci exp{—C 2 /T) + ai(0), where Ci and C 2 
are fit parameters. 


temperature phase diagram of the Random-Diluted TPM 
on the Bethe lattice must agree with the analytical solu¬ 
tion of the leaf-removal algorithm, that is also obtained 
on the random graph. From the leaf-removal analysis 
of Sec. |III| we already know that for all the concentra¬ 
tions of short-range plaquettes ag < 1 there is always 
a value > 0 such that for concentrations of long- 

range plaquettes q;l < a™'" the system is liquid at all 
tem pera tures. This happens because, as we discussed in 
Sec. |Tn] (see also Fig. [^, the only value of ag such that 
the related XOR-SAT problem becomes UNSAT (glass 
phase) for every > 0 is just ag = 1. Let us note that 
the phase diagram of our model in the plane (aL,T), 
i.e. for a fixed concentration ag of short-range plaque¬ 
ttes (see Fig. [^, has a remarkable similarity with the 
phase diagram in the plane (e, T) of^^, where e is an 
external field coupling different replicas of a TPM. Ac¬ 
cording tcJi^the critical line T'ii'(e) approaches the origin 
with infinite slop^i^, Tk ~ [log(e)]“^: we are going to 
show that the same happens in our case to when 

the concentration a^ of long-range plaquettes is sent to 
zero. We want to emphasize that plays a role analo¬ 
gous to e. The argument for the behaviour of Tx^OiL) in 
our system is rather simple and is exact on the random 
graph. FrorrP3 we know that in the UNSAT phase the 
ground states of a TPM on random graph have extensive 
energy. The energy of the system is by definition the 
sum over plaquette energies, so that when a ground state 
has extensive energy it means that there is an extensive 
amount of excited plaquettes in it. To have ground states 
with extensive energy is equivalent to have a minimum 
value, e““, for the energy per plaquette. In the TPM the 
concentration of defects at low temperatures behaves as 


c ^ and we can assume that in presence of a small 
amount of extra plaquettes this dependence is roughly 
the same, say c ^ with a ^ 1. Since the existence 

of a minimum value is equivalent to a minimum con¬ 

centration of excited plaquettes, also a minimum c™™ is 
fixed. From the constraint of a minimum allowed con¬ 
centration of excited plaquettes we can define a critical 
temperature Tk as Clearly the mini¬ 

mum value of the energy per plaquette and the minimum 
concentration of excited plaquettes must be proportional 
c™™ e™“. By looking at the behaviour of the energy in 

the ground states as a function of a presented iiP^l, it is 
reasonable to assume also for our Random-Diluted TPM 
on random graph that we have ~ a\, with 6^1. 
Putting together all the information collected above we 
can conclude that a^ = which in turn implies 


Tk 


a 2 
b log(aL) ■ 


(27) 


Eq. (271 allows us a good fit of the data in Fig. and 
accounts for the infinite slope of the curve T'if(aL) ap¬ 
proaching az, = 0. This infinite slope is signaling that 
around a^ = 0 the Random-Diluted TPM is sensitive to 
arbitrarily small perturbations, which induce the forma¬ 
tion of glass-phase. The same happens for the original 
two-dimensional Triangular Plaquette Model in presence 
of an external field e coupling replicaJi^. 


V. CONCLUSIONS 

In this paper we studied the thermodynamic proper¬ 
ties of the Triangular Plaquette Model in presence of ad¬ 
ditional plaquettes, namely we looked to what happens 
when a = M/N > 1, where M is the number of pla¬ 
quettes and N the number of spins. We have demon¬ 
strated that in the small-word lattice obtained by adding 
long-range plaquettes to the TPM the high-temperature 
expansion of the free-energy, which can be computed in 
the annealed approximation, has an entropic crisis. In 
the same model we find the numerical evidence of a first- 
order transition to a an ordered phase at T^, with a 
remarkably stable supercooled liquid phase at lower tem¬ 
peratures. The same phenomenology is found also when 
the additional interactions are short-range, although in 
this case there are corrections to the high-temperature 
expansion which very likely prevent the entropy crisis. 
Since the presence of an ideal glass transition is more 
likely when the entropy crisis takes place, we studied in 
more detail in the rest of the paper the model with long- 
range additional plaquettes, i.e. the Random-Diluted 
TPM. Our deepest investigation of the thermodynamic 
properties of the Random-Diluted TPM was carried on 
by means of the leaf-removal algorithm, usually applied 
to constraint satisfaction problem^^. The advantage of 
leaf-removal is that it allows to infer the thermodynamic 
properties of the model just analyzing the interaction net¬ 
work. The drawback is that the correspondence between 
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the properties of the interaction network and the ther¬ 
modynamics is exact only for random geometries. By 
means of the leaf-removal algorithm we obtained a tenta¬ 
tive phase diagram in the space of parameters ( 05 , 0 ;^), 
where Og is the concentration of plaquettes in the 2D 
triangular lattice while ol is the concentration of “long- 
range” plaqnettes. Such a phase diagram suggests that 
among the 2D plaquette models with different dilutions 
as the original TPlvfiS] (Os = 1 ) is the only one where 
the addition of any concentration of long-range plaque¬ 
ttes makes the thermodynamic non-trivial. That is why 
we say that the TPM is stochastically unstable: arbitrar¬ 
ily small perturbations of the Hamiltonian have dramatic 
effects on the thermodynamic properties of the model. 
Moreover, our results suggest that even in finite dimen¬ 
sions the parameter that controls this stochastic unstabil¬ 
ity is the ratio a between the number of plaquettes and 
the number of spins. The se con si derati ons are also com¬ 
patible with the results oP^Eol, jjJHEol jg jjow the 

TPM supports both dynamic and thermodynamic phase 
transition under the influence of arbitrarily small exter¬ 
nal fields. We find remarkable the similarity between 
the behaviour of the glass transition temperature Tk as 
function of Ol in our Random-Diluted TPM (on the ran¬ 
dom graph) and as function of e irP^ in both situations 
an infinitesimal amount of perturbation lifts the critical 
temperature Tk to finite values. We can conclude saying 
that further investigations on the stochastic stability of 
the TPM are mandatory. 


Let us recall that the correlation function in Eq. (Al) is 
the multispin correlation function (cti ... cr 3 m)TPM! where 
the 3m spins are randomly chosen with uniform probabil¬ 
ity on the 2D lattice. The multispin correlation function 
(cTi... o' 3 m)TPM is different from zero only when is pos¬ 
sible to find plaquettes in the TPM such that a set of 
one or more hyperloops can be formed which include all 
the 3m spins. Since in the TPM we know that both the 
magnetization (cri).pp„ = 0 and the two spin correlation 
function ((Ticrj)TPM = 0 zero, each of the hyperloops 
contains necessarily at least three spins. At the same 
time, since the 3m spins are chosen with uniform prob¬ 
ability on the lattice, the typical distance between any 
two of them is 0{\/N). This in turn means that the typ¬ 
ical the distance between any two spins which belong to 
the same hyperloop is also 0{'/N). Then, due to the 
fact that we consider the thermodynamic limit at fixed 
m, any hyperloop which contains the 3m spins of the 
the additional plaquettes contains also an infinite number 
of plaquettes of the TPM model. More precisely, since 
the hyperloops connects spins at a distance 0{VN), the 
number of plaquettes of the TPM in the hyperloop is 
0 (^dH/ 2 ), dn is the fractal dimension of the Sier- 

pinski gasket in D=2. One then has that the value of 
the correlation function between 3m spins chosen with 
uniform probability on the lattice is dominated at large 
N by the weight of the hyperloop which connects all the 
spins: (n ... Tm)TPM ~ [tanh(/3)]'^”^“^^. We can therefore 
conclude by noticing that for every m one has 
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lim 

N—>oo 



(n . . . Tm)^ 


A^™[tanh(/3)]'^"”'' =0. 

(A2) 


Appendix A: High-temperature expansion with 
additional plaquettes 

1. Random choice of spins in the new plaquettes: 
Random-Diluted TPM 

a. Triviality of the high-temperature expansion 

In order to say that the high-temperature expansion of 
Eq. [^in Sec. IIA is trivial we need to show that 


lim ( =0 

N—>-oo \ Ul 


V m 


(Al) 


b. Estimate of 


Making use of the argument discussed in App. (Alai, 


we show here that the partition function of the modified 
TPM of Sec. EH is self-averaging, namely that we have 
Z = Let us first write the expression of the high- 
temperature series of in a convenient way and then 
take the average over the disorder. We have that, for a 
given instance of the disorder, Z^ reads as: 
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, —/SffTPMlCT] —/3 J/tPm[s]+/3 Z]t-=^1 ^i(r)0'j(r)0'k(r)+Si(r)Sj(r)Sk(r) 


— 2 [cOSh(/3)] ^ ( J^(l “t” “t” '^i(r) (r) ^fc(r)] “t” tan}l(/?) (Jj ) 


where we have introduced the average: 


)ct,s — - 2 ^tpm 


E' 


, —/3(ffTPM[<7-]+f/’TPM[s]) 


(A3) 

(A4) 


It is then useful for what follows to define also: 


— ^TPM E ' 

cr 

= ^TPM E ' 


, —/3ffTPM[o’] 


—^Htpm[s] 


(A5) 


Now, in order to lighten the notation, is worthing to use the plaquette varia bles Tr = cFi(r)(Jj{r)<^k(r) and tr = 
Si{r)Sj(r)Sk(r) and use the symbol gr to represent the polynomial expression in Eq.(A3) corresponding to the plaquette 
“r”: 


Qr = tanh(/l)[rr + tr] + tanh(/3)^rrcrr 
so that we can write the high-temperature series of reads as 


22^[cosh(/3)]^^^+2^ A -f ^(5r) + '^{gvQp) + .. .y 

\ T r,p / 


(A6) 


(A7) 


the spins in each of the random plaquettes, is particularly simple and yields 


The average of the expression in Eq. (A7) over the disorder, which is represented by all the possible ways to choose 

(AS) 


Z2 = 22^[cosh(/3)]^^^ ^ 

so that, since in the limit A^ —)■ oo at fixed to we can write (^) AT™ = (ai-ZV)’", we have 


2"^" [cosh(/3)]"“^ < lim Z2 < 2"'''[coi 

N^OO 


/ oclN \ 

sh(^)]2^- • ■ • 9m).J 


The multiplaquette correlations in Eq. (A91 reads in turn 

/ m—k 




(i?!---!!™)=E u (n n 


J ‘3 


(A9) 


(AlO) 




In Eq.( AIO) the lowest degree correlations are those obtained taking the index k = 0, namely are of the kind (ri... Tm)c 


or {ti .. .Tra-k)cr{tm-k+i ■ ■ ■tm)s- The terms on the right hand side of Eq. (AIO) where a sing le correlation function 
appears are vanishing in the thermodynamic limit due to the same argument of App. ( |A I a ). With the same kind 
of arguments one can show that even the terms (ti ... Tm)cr{ti ... tm)s decay to zero in the thermodynamic limit fast 
enough to compensate the combinatorial prefactors, yielding finally the desired result 


lim {aN)"^{gi ... g„ 

Af —>00 


= 0 . 


From Eq. (A9) and Eq. (All) it follows finally 

^ = 22^[cosh(/3)]2“^+2^ 


(All) 

(A12) 


2. Additional plaquettes on a regular sublattice: 
smallest hyperloop via simulated annealing 


high-temperature expansion of Eq. ([^. In order to dis- 


In Sec. EH of the paper we mentioned a simulated 
annealing method to find the smallest hyperloop in the 
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cuss this method let us first rewrite the partition function 
as 

Z = ^ exp /3 ^ a^^ aj^ 

{a} I P 

= [cosh(/3)] " E n (l + tanh(/3)cri^crj^crfcp) 

(Ti,...,<7N P 

= [cosh(/3)]^ E n E [tanh(/3)cripCrjpCrfcJ"% 

(Ti,...,(T n P Up—0,1 

(A13) 

where the index p runs over all the plaquettes of the sys¬ 
tem, both the plaquettes of the original TPM and the ad¬ 
ditional plaquettes of the auxiliary sublattice introduced 


Hamiltonian in Eq. (A16). One introduces an effective 


in Sec. IIB In the last row of Eq. (A13| is convenient to 


explicitly write the product over the spins 

Z = [cosh(^)]^ E { n Nnh(/3)]' 

CTi,...,criv {ra} ^ p 

jj-^(Epeai"p)mod2 


= [cosh(/3)]^ Y \ n [tanh(/3)] 

• 1^1 + 


(AM) 


where now the index p labels the plaquettes around a 
given spin, p G di. The sum over v ariab les 
appearing in the last rows of Eq. (A14| represents the 
sum over all possible collections of plaquettes, either 
forming or not an hyperloop. Within a certain collec¬ 
tion of plaquettes the one labeled with p is present when 
rip = 1 and absent when rip = 0. Let us stress that each 
collection of plaquettes which does not form an hyper¬ 
loop, i.e. an assignment of the variables rip such that, 
at least for one i, we have iJ2p^di ’^p)mod 2 = 1, does 
not contribute to the sum in Eq. ( |AI4[ ). That is why, in 
order to seek non-trivial terms of Z, we need to find hy¬ 
perloops. A hyperloop correspond therefore to a choice 
of {rip} such that for each spin we have 


E’ 

\.pGdi 


= 0 


(A15) 


mod2 


Therefore, in order to find hyperloops, one can look for 
the ground states of the dual model defined by the fol¬ 
lowing energy function: 


N 


Hdual(iV) = EE’ 

i—1 p^di 


(A16) 


The simulated annealing method comes at this stage as 
the most natural one to seek for the ground states of the 


inverse temperature parameter j3eS = T~f[ and then sam¬ 
ples configurations according to the Boltzmann measure 
exp(—/3eff'Hduai(A^))j while slowly decreasing T^s, until a 
configuration with if = 0 is found. In order to find the 
smallest hyperloop we have realized this simulated an¬ 
nealing search varying the size A^ of a TPM with an aux¬ 
iliary sublattice made of side 3 triangles and with open 
boundary conditions. Varying N we considered always 
lattices with a side that was a multiple of the incommen¬ 
surate sublattice cell side, namely we considered triangu¬ 
lar lattices of rhomboidal shape and side L — 6,9, 12,.... 
The result of our study is that the smallest hyperloop ap¬ 
pears in lattice of side L = 12 made of plaquettes of 
the TPM and 4^ plaquettes of the incommensurate sub¬ 
lattice. Such a hyperloop, which is represented in Fig[^ 
is made by 54 plaquettes of the TPM model and 10 pla¬ 
quettes of the incommensurate sublattice. Let us notice 
that, as argued above in App. (Ala), the hyperloop we 


have found is symmetric with respect of the three sym¬ 
metry axes of the lattice. 


Appendix B: Leaf-removal exact solution for the 
Random-Diluted TPM on random graph 

For a random regular graph the action of the “leaf- 
removal” algorithm is represented in the thermodynamic 
limit by an infinite set of differential equationJ^ for the 
connectivities of spins, ni(t), where t is the reduced time 
t = n/N, with N the total number of spins in the system 
and n the number of iterations of the leaf-removal algo¬ 
rithm. Because the maximum possible value taken by n 
is N the reduced time is in the interval [0,1]. Our Ran¬ 
dom Diluted-TPM differs from the XOR-SAT because 
for the latter the initial connectivity ni{0) is Poissonian, 
while in our case is the mixed Poissonian/binomial dis¬ 
tribution of Eq. (23) in Sec. |III| o f the main text. The 
probability distribution in Eq. ( |23[ ) is the only ingredient 
of the original Random Diluted-TPM left when the model 
is studied on the random regular graph. While the be¬ 
haviour of ni{t; as, ol) can be studied analytically on the 
random graph, for the original RD-TPM model this can 
be studied only numerically. The differential equations 
for the evolution of the connectivity under leaf-removal 
are 


niit) — —^£i+<Jfo+ 


3(q:s + ttL -t) 


[{£ + - £ni(t)], 

(Bl) 


where b(t) = (I — t/{as -I- The general solution 

for £ > 2 is: 


ntit) = b^\t)Yne+kiO) 

p=0 


[-b^{t)Y 


1 




r=0 


"! (k — r)\ 


(B2) 

By assuming the distribution in Eq. (23) as initial condi¬ 
tion n^+fc(0) and plugging it into Eq. (B2) one finds for 
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> 2 : 


ni{t) = 


[3aLb'^{t)Y e 


l „-3ai,b=(t) 3 


e. 


N 


Y,Cr,iBr{g{t)), (B3) 


r =0 


where Bj. (x) is the Bell polynomial of order r in the vari¬ 
able X, g{t) = 3q;l(1 — 6 ^(t)) and the coefficients Cr^e 
are: 


(3^ ^ 


P3 


(3aL)= 


Ci,f = /^ + (2€ - 1) 


P2 


C 24 = 


C^,i = 


3q!l 
P2 


{3aLY 

P3 

(3aL)3- 


3(^-1) 


(3aL)2 

P3 


+ (3^2 -6£ + 2) 


P3 


(3aL)' 


(3aL)' 


The probability that a certain spin is a leaf at the itera¬ 
tion t, ni(t), is then: 


ni{t) = 


ni( 0 ) + ds 

Jo 


4n2(s) 


3(as -I- CKL - s) 


- 1 


1 


62(s) 

(B5) 

When the probability to find a leaf vanishes before all the 
spins are eliminated, i.e. when ni(t) = 0 with t < 1 , the 
system has a finite core. On the contrary when ni (t) = 0 
only at t = 0 there is no core. Studying the behaviour of 
ni(t) at different values of as and ul it is then possible 
to locate in the parameter space the line (as”®,a™''®). 
In order to find the SAT/UNSAT transition one needs 
then to calculate the number of spins and plaquettes in 
the core, when it is present. The concentration of spins 
in the core reads 


i=2 


(B 6 ) 


and by plugging the definition of niit) from Eq. (B2) into 
Eq. (B 6 ) one gets: 

3 


ndt) = 


o-g{t) 

~w~ 


-ICo{t) -f ^ K.r(t){-g{t) -f e^^*'>Br-{g{t))) 


r—0 


(B7) 

where we have defined g{t) = 3aLb"^{t) and the coefh- 
cients K-rif) read as: 


Pi 

3aL 


ICoit) = po + f{t) 

+ {fit) + fit)) 


P2 


(3aL)' 

P2 


+ 2 


P3 


(3aL)3j 

P3 


(3aL)2 ^3ai)3j 

{fit) + V\t) + fit))^^ 


lCi{t) = ^- 

Ml 

+2/(t) 
+ {ft) 


P2 


+ 2 


P3 


(3aL)2 
P2 

(3ai)2 

£fit))3 


- 3 


(3ai)3 
P3 


(3aL)3j 

P3 


(3aL)3 


JC2it) = 
left) = 


P2 


(3aL)2 

P3 

(3aL)3 


-3 


P3 


(3aL)3 




(3aL)3 


(B 8 ) 


The number of plaquettes left in the system at time t 
in the leaf-removal algorithm is then simply provided by 
rricit) = Us + — t. For each point (a^, af in the pa¬ 

rameter space one must look for the time t* at which the 
leaf-removal algorithm stops and then, if ndt*) > 0 , con¬ 
sider the ratio 7 = mdt*)/ndt*)'- the static transition 
line is identified by 7 = 1 . 


(B4) 


Appendix C: Leaf removal and partial traces 


Let us justify here some of the arguments in Sec. |IIIB| 
of the main text. We briefly explain why in the case 
when leaf-removal eliminates all the spins from the in¬ 
teraction network then also the partition function can be 
summed by means of partial traces without giving rise to 
any close diagram: this is the case when Z has the simple 
expression in Eq. (251 of Sec. Ill B of the main text. The 


method of partial traces consist in summing the partition 
function starting from the open boundaries of the lattice. 
The necessary condition for this method to work is that 
each step there is at least one spin belonging to a single 
plaquette. Let us explain how the iterative summation 
algorithm works: 

1. Look for a spin which appears in a single plaquette, 
say (To- 

2. Sum over values of the spin erg in the partition func¬ 
tion: 

/ M,+Mi \ 

Z = ^ exp ( P ^ j 

<70,...,<Tiy \ fj,= l / 


= cosh(,5) E [1 + no(^)Cii(^)cr 2 (^) tanh(/3)] • 

fTo,... 

( Ms+Ml 


(Cl) 


3. Check if either or cr 2 (,^) in Eq. (Cl), or both, 
are again participating to a single interaction 'y d v 
after the removal of v by summation. If for instance 
is appearing only in (Ti(j^)(Tj(.y)CTfc(.y), go back 
to point 2). If both cfk^u) and c! 2 {u) are participating 
to more than one interaction, go back to point 1 ). 
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4. If and only if the above recursion can be iterated 
summing over all the spin of the system the parti¬ 
tion function is the one of Eq. (251. 


From the description of the method of partial traces, 
it is clear that it is exactly the iterative scheme of leaf- 
removal. The situation when leaf-removal leaves a finite 
core are the situations where closed diagram arise in 2 , 
so that 2, cannot be simply summed by means of par¬ 
tial traces. It is worthing to recall that in order to use 
leaf-removal or partial traces one must choose the cor¬ 
rect boundary conditions: this is due to the determinis¬ 
tic nature of these decimation algorithms. Consider for 
instance the pure Triangular Plaquette Model with N 
spinal the partition function is exactly known to be 
Z = 2^[cosh(/3)]^. Nevertheless, in the case of 'peri¬ 
odic boundary conditions there are no leaves (no spins 
appearing in a single plaquette), so that both the leaf- 
removal algorithm and the partial traces iterations can¬ 
not be started: it looks like the whole system is a core 
of dimension N. On the contrary, if one considers open 
boundary conditions, one finds from partial traces that 
Z = 2'^[coshlog(/3)]'^, apart from corrections negligi¬ 
ble in the thermodynamic limit. Open boundary con¬ 
ditions are therefore the correct choice if one wants to 
use this kind of algorithms to study the thermodynam¬ 
ics of the system. Moreover, in a phase diagram like 
the one of Fig. open and periodic boundary conditions 
correspond to the same points, because they differ for a 
sub-extensive number of plaquettes. 


Appendix D: Belief-propagation equations for the 
Random-Diluted TPM 


As mentioned in the main text in Sec |IV[ in order to 
represent the two different kind of plaquettes in the sys¬ 
tem, it is convenient to introduce two different cavity 
fields, and u^^i, respectively for the “short” and 

“long” range plaquettes. These fields allow one to write 
the marginal probability distribution of the spin i when 
all interactions around it but 7 are removed, with 7 rep¬ 
resenting respectively a short or long range plaquette: 


Pv{<^i) 2 cosh{l3vj^i) 2 cosh{f3uj^i) 

(DI) 

In order to write in a clear way the Belief Propagation 
(BP) equations, we need to introduce also the couple of 
cavity fields and Ui^-y, which represent the effec¬ 

tive field on Oi when only the plaquette 7 is removed, 
respectively when 7 is short and long range. The belief 


propagation equations for our model read then: 

riL — l ris 

pedj\'r Pedj\'r 


riL Tig —1 



P&dj\j 

/ \ 

tanh“^ j 

tanh(/3) tanh(/3Mj_>.y) 


\ jed'rA / 

/ \ 

tanh“^ j 

tanh(/3) tanh(/3{;j_>.^) j 

P ’ 

\ jed'rA / 


(D2) 


where til and rig denotes respectively the number of long 
and short range plaquettes attached to each spin j. Let 
enclose the belief propagation equations in Eq. (D2 1 in 
the expression 


^7—fz ‘A ({'Uy5_^J . 


(D3) 


The population dynamics algorithm is realized starting 
with a sufficiently large sample of values for each of the 
two fields u and v, randomly initialized with flat distribu¬ 
tion in the interval [—1, 1 ]. A random sequential update 
of the values in the two arrays is realize according to 
the BP equations in Eq. (D2). The numbers of long, 


and short range plaquettes, rig, attached to each spin and 
necessary for each iteration step of the algorithm are ran¬ 
dom variables extracted according to the distributions of 
Eq. (22 1 in the text. As is clear from the first two lines 


of Eq. (D2), at each iteration step one also needs the 


excess degree distributions Pexc and Pexc defined as fol¬ 
lows. If we already know that the long-range plaquette 
7 is attached to the spin cji, pexc is the probability that 
ul — 1 other long-range plaquettes are attached to cr^. 
The same is true for the definition of pexc when we know 
that 7 is a short range plaquette. Such distributions read 
respectively 


PexcinL - 1) = 


riLpiriL) 

{nr) 


Pexcij^s 1) — 


spjng) 

{ng) 


(D4) 

On the Bethe lattice our model has both a dynamic phase 
transition, at T^, and thermodynamic transition, at Tk- 
While the dynamic ergodicity breaking at Td disappears 
in interaction networks with with finite loops, the ideal 
glass transition at Tk may survive. The dynamical tran¬ 
sition temperature Td correspond to the formation of an 
exponentially large number of metastable states sepa¬ 
rated by extensive barriers and such that the system, 
when initialized in one of this states is trapped within 
it. In term of the cavity equations this phenomenon can 
be recognized by introducing a couple of auxiliary fields 
Uo-=±i (an Vcr=±i) for each type of cavity field {u and 
v). The cavity fields represent the value of the field 
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on a conditioned to the knowledge of the value taken by 
this spin, either a = 1 or a = —1. If the populations 
of Ucr=i and Ucr=-i, whic h ar e calculated according to 
the equations below, Eq. ( |D5[ ), at equilibrium are such 
that the V{u) = Vi{ui) = where V, Vi and 

V-i are the probability distributions respectively of u, 
M_i and u_|_i, it means that the system is in the simple 


paramagnetic state. On the contrary when the two distri¬ 
butions 7^1 (ui) and V-i{u-i) become different from the 
the distribution V{u) of the equilibrium field, it means 
that the system “remembers” the initial condition and 
the effective field around a certain spin favors the values 
taken by such spin in the initial condition: ergodicity is 
dynamically broken. The equation to recursively update 
the distributions Va(ua) is the following. 


'Pa{Ua\u) = ^ J7Pexc("lL(i))p(ns(i)) J 


{mL{i),ns{i)} *=1 


2 mL(i)ns{i) 

nn n du^ dv'^V {u^)V ) 

i—1 j k 




E 






nn n 


2 na(i) 

i—lj k 


2cosh(/3uf) 2cosh(,5r;'=) 

S{u^ - 


(D5) 


where ns{i) is drawn from in Eq. ( |22[ ) while rriLii) is 
drawn from pea:c in Eq. ( |D4[ ) . The same kind of equations 
holds for Vcriva-lv), just with 71 ^( 1 ) (drawn from in 
Eq. (l22|_in place of rnL{i) and ms{i) (drawn from 
in Eq. ( |D4[ ) in place of ns{i). 

The iteration step of the population dynamics accord¬ 
ing to Eq. (D5) proceed as follows: 


1. Choose an element to update in the population of 
fields Mcr, which is equivalent to say: choose ran¬ 
domly a spin cr in the lattice. It is assumed that 
we are interested in the cavity field on a that is ob¬ 
tained by removing all the plaquettes but one, say 
the plaquette 7 . When studying the distribution of 
u„ we know that 7 is a long range plaquette. 


2. Consider the spins Ui which are interacting with a 
through the plaquette 7 . Extract then the number 
of long-range, mL{i) — 1 , and short-range, ns{i), 


plaquettes attached to each of the spins <Ji. 


3. Compute the cavit y field u according to the func¬ 
tion in Eq. (D3| from the cavity fields {u^ 


4. Choose the values of spins Oi according to 
equilibrium measure for a given value of 
(7, namely with probability p((7i, cr 2 |cr) = 

y-r2 y-rmi.(7 -r-rn^(7 

Z({u^,v^}) life 2cosh(^'U-?) 2 cosh(/3i!^) ■ 


5. 


For each cavity held (v^) consider the attached 


6 . From the set of {m^. update the held Ucr (u^) 
according to J-. 


Once the population dynamic algorithm is converged 
and the stationary probability distributions also for the 
conditioned cavity helds Mg- and Va- is know, one can cal¬ 
culate from it the free-energy within a single metastable 
state according to the following formula. 


/meta = aL{f l{u,V ,Uc, ,Va)) + as{fs{u,V ,Ua ,Vc,)) - ^ Qn{cts,aL){n - l){f^"-\u,V,Ua,Va)), 

n 


(D 6 ) 


where we have called here Qn{cis,OiL) the mixed Poisson/binomial probability distribution of spin connectivity 
dehned in Eq. (23) as ni{as, ocl)- The free-energy per plaquette reads, in the case of a long range plaquette, as 
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{fL{u,V,Ua,Va)) = ^ ^ f[pexc{mL{i))pins{i)) f 

{(TiljiGOA {m L {i) ,n s (i)} ^—1 


2 mL{i)ns{i) 

nn n du^dv^V{u^)V{v'‘) 


i—1 j k 




cr- 3 mL{i)n^(i) gj 

A£!z^ff TT fr ^^_^_ 

Z{{ui,v^}) 2cosh(/3M'5) 2cosh(/3i;^) 

2 rriLii) ns{i) 

n n duiTMMn n 

i—1 j k 




(D7) 


The expression in Eq. (D7) that has to be consistently modified for a short range plaquette. The free-energy per spin 
is then 


nn du^dv^V{u^)V{v^)) 

j k 


{^} 

mL{i) ris{i) 


^ ruL 


nn 


Z({uf,w^}) -fy*- Y 2cosh(/3M*) 2cosh(/3w'=) 

J — 1 k —1 

logZ<,(K,z;^}). 


(D8) 


The two partition functions Z^ and Z^ appearing respectively in Eq. (D7) and Eq. (D81 are define as 


ZA{{ui.,v’;^}) = Y^ 

{^i} 


3 mi(i)"=(i) 

/3nLi'^. ■TT TT T-r e e 


n n rr 2cosh(/3Mf)2cosh(/3'y^) 


z.({<,<})=^nn 


a ^j3v^<7 


cr 2cosh(/3uf) 2cosh(/3n^) 

The free energy in the paramagnetic phase is the same as in the high-temperature expansion, 

/para = -^"^log(2) + ctl) logcosh(/3). 

The configurational entropy can be finally obtained as: 

^ Pi^fmeta. fpa.ra) 


(D9) 

(DIO) 

(Dll) 
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